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Can a gas of spin-up and spin-down fermions become ferromagnetic due to 
repulsive interactions? This question which has not yet found a definitive 
theoretical answer was addressed in an experiment with an ultracold two- 
component Fermi gas. The observation of non-monotonic behavior of lifetime, 
kinetic energy, and size for increasing repulsive interactions provides strong 
evidence for a phase transition to a ferromagnetic state. It implies that itin- 
erant ferromagnetism of delocalized fermions is possible without lattice and 
band structure and validates the most basic model for ferromagnetism intro- 
duced by Stoner. 

Magnetism is a macroscopic phenomenon with its origin deeply rooted in quantum me- 
chanics. In condensed matter physics, there are two paradigms for magnetism: localized spins 
interacting via tunnelling, and delocalized spins interacting via an exchange energy. The latter 
gives rise to itinerant ferromagnetism which is responsible for the properties of transition metals 
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like cobalt, iron and nickel. Both kinds of magnetism involve strong correlations and/or strong 
interactions and are not yet completely understood. For localized spins, major open questions 
include the interplay of magnetism with d-wave superfluidity and frustrated spin materials. For 
itinerant ferromagnetism {1-7), phase transition theories are still qualitative. 

We implement the Stoner model, a textbook Hamiltonian for itinerant ferromagnetism (<S), 
using a two-component gas of free fermions with short-range repulsive interactions, which can 
capture the essence of the screened Coulomb interaction in electron gases {8). However, there 
is no proof so far that this simple model for ferromagnetism is consistent when the strong 
interactions are treated beyond mean-field approaches. It is known that this model fails in one 
dimension where the ground state is singlet for arbitrary interactions, or for two particles in 
any dimension (5). Here, cold atoms are used to perform a quantum simulation of this model 
Hamiltonian in 3D and show experimentally that it leads to a ferromagnetic phase transition (2). 
This model is also realized in helium-3 (9), but it turns into solid and not into a ferromagnetic 
phase at high pressure. It has also been applied to neutrons in neutron stars {10). 

So far, magnetism in ultracold gases has been studied only for spinor {11, 12) and dipo- 
lar {13) Bose-Einstein condensates. In these cases, magnetism is driven by weak spin-dependent 
interactions which nevertheless determine the structure of the condensate due to a bosonic en- 
hancement factor. In contrast, here we simulate quantum magnetism in a strongly interacting 
Fermi gas. 

An important recent development in cold atom science has been the realization of super- 
fluidity and the BEC-BCS crossover in strongly interacting two-component Fermi gases near a 
Feshbach resonance {14). These phenomena occur for attractive interactions for negative scat- 
tering length and for bound molecules (corresponding to a positive scattering length for two 
unpaired atoms). Very little attention has been given to the region for atoms with strongly 
repulsive interactions. One reason is that this region is an excited branch, which is unstable 
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against near-resonant three-body recombination into weakly-bound molecules. Nevertheless, 
many theoretical papers have proposed a two-component Fermi gas near a Feshbach resonance 
as a model system for itinerant ferromagnetism (15-22) assuming that the decay into molecules 
can be sufficiently suppressed. Another open question which has not been addressed is the 
possibility of a fundamental limit for repulsive interactions. Such a limit due to unitarity or 
many-body physics may be lower than the value required for the transition to a ferromagnetic 
state. We show that this is not the case, and that there is a window of metastability where the 
onset of ferromagnetism can be observed. 

A simple mean-field model captures many qualitative features of the expected phase transi- 
tion, but is not adequate for a quantitative description of the strongly interacting regime. The 
total energy of a two-component Fermi gas of average density n (per spin component) in a vol- 
ume V is given by EF2Vn[^{{l + r]f/^ + (1 - r]f'^} + |^^Fa(l + ^) (1 - r])] where Ep is the 
Fermi energy of a gas, kp the Fermi wavevector of a gas, a the scattering length characterizing 
short-range interactions between the two components, and r/ = ^n/n = (ni — 712) / (ni + 77,2) 
magnetization of the Fermi gas. The local magnetization of the Fermi gas is non-zero when the 
gas separates into two volumes, where the densities rii and 712 of the two spin states differ by 
2 A?^. Note that we study an ensemble in which the number of atoms in each spin state is con- 
served. This is equivalent to a free electron gas at zero external magnetic field where the total 
magnetization is zero. The interaction term represents any short-range spin-independent poten- 
tial. When the gas is fully polarized, it avoids the repulsive interaction, but increases its kinetic 
energy by a factor of 2^/^. The phase transition occurs when the minimum in energy is at non- 
zero magnetization (Fig. lA) at kpa — 7r/2. This onset was discussed in the context of phase- 
separation in a two-component Fermi gas (15-18). Fig. IB shows several consequences of the 
phase transition for a system at constant pressure. First, for increasing repulsive interactions, the 
gas expands, lowering its density and Fermi energy; kinetic energy is therefore reduced. When 
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the gas enters the ferromagnetic phase, kinetic energy increases rapidly due to the larger local 
density per spin state. Furthermore, the volume has a maximum value at the phase transition. 
This can be understood by noting that pressure in our model is (2/3)£^jti„/\^ + Eint/V, where 
Ekin is kinetic energy and Ei^t interaction energy. At the phase transition, the system increases 
its kinetic energy and reduces its interaction energy, thus reducing the pressure. This maximum 
in pressure at constant volume turns into a maximum in volume for a system held at constant 
pressure, or in a trapping potential. We have observed three predictions of this model: the onset 
of local magnetization through the suppression in inelastic collisions, the minimum in kinetic 
energy, and the maximum in the size of the cloud. These qualitative features are generic for the 
ferromagnetic phase transition and should be present also in more advanced models (19). 

We start with an atom cloud consisting of an equal mixture of ^Li atoms in the lowest 
two hyperfine states, held at 590 G in an optical dipole trap with additional magnetic confine- 
ment (25). The number of atoms per spin state ~ 6.5 x 10^ corresponds to a Fermi temper- 
ature Tp of ~ 1.4 fxK. The effective temperature T could be varied between T/Tp = 0.1 and 
T/Tp = 0.6 and was determined right after the field ramp by fitting the spatial distribution of 
the cloud with a finite temperature Thomas-Fermi profile. Note that kp is the Fermi wavevector 
of the non-interacting gas calculated at the trap center. Applying the procedure discussed in 
Ref. (24) to repulsive interactions, we estimate that the real temperature is ~20% larger than 
the effective one. The effective temperature did not depend on kpa for kpa < 6. At higher 
temperatures, additional shot-to-shot noise was caused by large fluctuations in the atom num- 
ber. From the starting point at 590 G, the magnetic field was increased towards the Feshbach 
resonance at 834 G, thus providing adjustable repulsive interactions. Due to the limited lifetime 
of the strongly interacting gas, it was necessary to apply the fastest possible field ramp, limited 
to 4.5 ms by eddy currents. The ramp time is approximately equal to the inverse of the axial 
trap frequency (25) and therefore only marginally adiabatic. Depending on the magnetic field 
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during observation, either atoms or atoms and molecules were detected by absorption imaging 
as described in Fig. SI (25). 

The emergence of local spin polarization can be observed by the suppression of (either 
elastic or inelastic) collisions, as the Pauli exclusion principle forbids collisions in a fully po- 
larized cloud. We monitor inelastic three-body collisions which convert atoms into molecules. 
The rate (per atom) is proportional to f{a,T)nin2 or /(a,T)n^(l — rj'^) and is therefore a 
measure of the magnetization rj. For kpa <S 1, the rate coefficient /(a, T) is proportional to 
max(T, Tp) {26). This rate can be observed by monitoring the initial drop in the number of 
atoms during the first 2 ms after the field ramp. We avoided longer observation times since the 
increasing molecule fraction could modify the properties of the sample. 

Fig. 2 shows a sharp peak in the atom loss rate around k°pa ~ 2.5 at T/Tp = 0.12 indicating 
a transition in the sample to a state with local magnetization. The gradual decrease is consistent 
with the inhomogeneous density of the cloud where the transition occurs first in the center. The 
large suppression of the loss rate indicates a large local magnetization of the cloud. 

The kinetic energy of the cloud was determined by suddenly switching off the optical trap 
and the Feshbach fields right after the field ramp and then imaging state 1 1) atoms at zero field 
using the cycling transition after a ballistic expansion time of A^^j = 4.6 ms. The kinetic energy 
was obtained from the Gaussian radial width Ux as Ekin = where m is the mass of the 

^Li atom. Fig. 3 demonstrates a minimum of the kinetic energy at k°pa ~ 2.2 for the coldest 
temperature, T/Tp = 0.12, nearly coinciding with the onset of local polarization. The peak in 
the atom loss rate occurs slightly later than the minimum of kinetic energy, probably because 
/(a, T) increases with a (22). Since the temperature did not change around kpa ~ 2.2, the 
increase in kinetic energy is not caused by heating, but by a sudden change in the properties of 
the gas, consistent with the onset of ferromagnetism. The observed increase in kinetic energy is 
~20 % at T/Tp = 0.12, smaller than the value (2^/^ - 1) = 0.59 predicted for a fully polarized 
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gas. This discrepancy could be due to the absence of polarization or partial polarization in 
the wings of the cloud. Also, it is possible that the measured kinetic energy of the strongly 
interacting gas before the phase transition includes some interaction energy if the Feshbach 
fields are not suddenly switched off. For the current switch-off time of ~ 100 /is, this should be 
only a 5% effect, but the magnetic field decay may be slower due to eddy currents. 

Finally Fig. 4 shows our observation of a maximum cloud size at the phase transition, in 
agreement with the prediction of the model. The cloud size may not have fully equilibrated 
since our ramp time was only marginally adiabatic, but this alone cannot explain the observed 
maximum. 

The suppression of the atom loss rate, the minimum in kinetic energy, and the maximum 
in cloud size show a strong temperature dependence between T/Tp of 0.12 and 0.22. As the 
properties of a normal Fermi gas approaching the unitarity limit with k^a » 1 should be 
insensitive to temperature variations in this range, this provides further evidence for a transition 
to a new phase. 

At higher temperature (e.g. T/Tp = 0.39 in Fig. 3), the observed non-monotonic behavior 
becomes less pronounced and shifts to larger values of k^a for 3 < k^a < 6. For all three 
observed properties (Figs. 2-4), a nonmonotonic behavior is no longer observed at T/Tp = 
0.55 (27). One interpretation is that at this temperature and above, there is no phase transition 
any more. Note that in a mean-field approximation, a ferromagnetic phase would appear at all 
temperatures, but for increasing values of k^a. Our observations may imply that the interaction 
energy saturates around k^a 5. 

The spin-polarized ferromagnetic state should not suffer from inelastic collisions. However, 
typical lifetime were 10 - 20 ms, probably related to a small domain size (see below) and three- 
body recombination at domain walls. 

We were unsuccessful in imaging ferromagnetic domains using differential in-situ phase- 
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contrast imaging (28). A noise level of 5'/A'^ ~10 suggests that there were at least 100 domains 
in a volume given by our spatial resolution of ~ 3 fim and the radial size of the cloud. This 
implies that the maximum volume of the spin domains is ~ 5 /im^, containing ~ 50 spin- 
polarized atoms. We suspect that the short lifetime prevented the domains from growing to 
a larger size, and eventually adopting the equilibrium texture of the ground state, which has 
been predicted to have the spins pointing radially outward, like a hedgehog (20, 22). All our 
measurements are sensitive only to local spin polarization, independent of domain structure and 
texture. 

The only difference between our experiment and the ideal Stoner model is a molecular 
admixture of 25 % (Fig. 4). The molecular fraction was constant for k^a > 1.8 for all tem- 
peratures and therefore cannot be responsible for the sudden change of behavior of the gas at 
kpa ~ 2.2 for the coldest temperature T/Tp — 0.12 . This was confirmed by repeating the 
kinetic energy measurements with a molecular admixture of 60 %. The minimum in the kinetic 
energy occurred at the same kpa within experimental accuracy. 

Before we can compare the observed phase transition at k^a ~ 2.2 to the theoretical pre- 
dictions, we have to replace the ideal gas kp by the value for the interacting gas, which is 
smaller by ~ 15% because of the expansion of the cloud (Fig. 4), and obtain a critical value 
for kpa ~ 1.9 ± 0.2. At T/Tp — 0.12, the finite temperature correction in the critical value 
for kpa is predicted to be less than 5% (19). The observed value for kpa is larger than the 
mean-field prediction ofn/2 and the second order prediction of 1.054 at zero temperature (19). 
Depending on the theoretical approach, the phase transition has been predicted to be first or sec- 
ond order. This could not been discerned in our experiment due to the inhomogeneous density 
of the cloud. 

Ref. (19) speculated that earlier experiments on the measurement of the interaction en- 
ergy (29) and RF spectroscopy of Fermi gases (30) showed evidence for the transition to a 
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ferromagnetic state at or below kpa — 1. This interpretation is ruled out by our experiment. 

Our work demonstrates a remarkable asymmetry between positive and negative scattering 
length. Early work {15) predicted that for /c^lal = 7r/2, both an attractive and a repulsive Fermi 
gas become mechanically unstable (against collapse, and phase separation, respectively). In an 
attractive Fermi gas, however, the mechanical instability does not occur (due to pairing (31)), 
in contrast to our observations in a repulsive Fermi gas. This suggests that the maximum total 
repulsive energy (in units of ?>/b{2Vn)EF) is larger than the maximum attractive energy of 
0.59 {32) realized for infinite a {23). 

Heisenberg's explanation for ferromagnetism was based on exchange energy, i.e. the Pauli 
principle and spin-independent repulsive interactions between the electrons. However, it re- 
mained an open question, what other "ingredients" were needed for itinerant ferromagnetism. It 
was only in 1995 {6, 7), that a rigorous proof was given that, in certain lattices, spin-independent 
Coulomb interactions can give rise to ferromagnetism in itinerant electron systems. Our find- 
ing implies that Heisenberg's idea does not require a lattice and band structure, but applies 
already to a free gas with short-range interactions. Our experiment can be regarded as quantum 
simulation of a Hamiltonian for which even the existence of a phase transition was unproven. 
This underlines the potential of cold atom experiments as quantum simulators for many-body 
physics. 
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Figure 1 : Ferromagnetic phase transition at T=0, according to the mean-field model described 
in the text. The onset of itinerant ferromagnetism occurs when the energy as a function of 
magnetization flips from a U-shape to a W-shape (A). Figure (B) shows the enthalpy, volume 
and kinetic energy (normalized to their values for the ideal Fermi gas), and magnetization as a 
function of the interaction parameter kpa. Note that kpis defined by the density of the gas. The 
dotted line marks the phase transition. 
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Figure 2: Atom loss rate as a probe for local spin polarization, for different temperatures, (a) 
T/Tf = 0.55 (dashed curve), (b) T/Tp = 0.22 (dotted curve), and T/Tp = 0.12 (solid black 
curve). The curves are guides to the eye, based on the assumption of a loss rate which saturates 
for increasing a in the normal state. The shaded area around the phase transition at T/Tp = 0.12 
highlights the same region as in Figs. 3 and 4. 
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Figure 3: Kinetic energy of a repulsively interacting Fermi gas determined for different inter- 
action parameters kpa and temperatures. The measured kinetic energy is normalized by the 
Fermi energy Ep of the noninteracting Fermi gas at T=0, calculated at the trap center with the 
same number of atoms per spin state. Each data point represents the average of three or four 
measurements. 
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Figure 4: Maximum in volume at the phase transition (A) Axial size and chemical potential of 
the cloud for various temperatures. The chemical potential /i is determined from the measured 
cloud size, cr^ as /i = \mujz'^az'^ (B) Number of particles including both atoms and molecules 
right after the field ramp. This result shows that 25 % of atoms were converted into molecules 
during the field ramp, and this fraction stayed constant for k°pa > 1.8, where the phase transition 
was reached. This molecule fraction was independent of temperature. 
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Supporting materials: 

Itinerant Ferromagnetism in a Fermi Gas of Ultracold Atoms 



Materials and Methods 

Preparation of the ultracold ^Li cloud The first step is the production of a spin-polarized 
Fermi gas in the \F = 3/2, mj? = 3/2) state by sympathetic cooling with bosonic ^^Na atoms in 
a magnetic trap as described in ref (SI). The ^Li cloud was then loaded into a deep optical dipole 
trap with a maxium power of 3W and radial trap frequency of ~3.0 kHz, followed by an RF 
transfer into the lowest hyperfine state |F = l/2,mir = l/2). Additional axial confinement was 
provided by magnetic fields. An equal mixture of |1) and |2) spin states (corresponding to the 
\F — 1/2, mp — 1/2) and |F = 1/2, mp — —1/2) states at low magnetic field) was prepared 
by a Landau-Zener RF sweep at a magnetic field of 590 G, followed by 1 s for decoherence and 
further evaporative cooling at 300 G. Finally, the optical trapping potential was adiabatically 
reduced over 600 ms, and the field increased back to 590 G. The trap had a depth of 7. 1 /jK and 
was nearly cigar shaped with frequencies — ^ 300 Hz and Uz ~ 70 Hz. 

Supporting online text 

Estimation of the maximum total repulsive energy Full phase separation at zero tempera- 
ture requires a total repulsive energy of (2^/^ — 1) = 0.59 in units of 3/5{2Vn)EF. At finite 
temperature T, one has to add TS where S ~ {2Vn)kBln2 is the entropy difference between 
the two phases. Our tentative observation of a ferromagnetic phase at T = 0.39Tp implies a 
repulsive energy of ~ 1.04 assuming full phase separation, larger than the maximum attraction 
energy of 0.59. 
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Fig. SI: (A) The schematic shows the time sequence of the experiment. The sample was 
exposed to the magnetic field of interest for - 14 ms and analyzed in-situ for loss measurement 
or after 4.6 ms time-of-flight for the measurement of kinetic energy and the axial size of the 
cloud. The Feshbach fields were suddenly switched off at a rate of lG//is, preventing the 
conversion of interaction energy into kinetic energy during the expansion. (B) This absorption 
image shows the |1) component of the cloud trapped at 812 G (left), and after 4.6 ms ballistic 
expansion imaged at zero field (right). The field of view is 840/umx 550/um. The magnetic field 
ramp was limited by eddy currents to 4.5 ms. Spectroscopic measurements of the magnetic field 
showed that the field was trailing behind the current which was controlled with a time constant 
faster than 1 ms. 
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